How large can be SDW and CDW amplitudes in underdoped cuprates 
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Self-consistent calculation of spin (charge) density wave order (SDW/CDW) parameters have 
been performed for bilayered cuprates on the basis of a singlet correlated band model. Evolution 
of the Fermi surface in the strongly underdoped regime is described by using a two-band approach. 
The smooth development of the pseudogap formation temperature is explained from underdoped to 
overdoped states and the Fourier amplitudes < > (spin) and < Cq > (charge) modulations have 
been calculated. We have found a maximum of the incommensurability for doping 0.09 -7- 0.11 holes 
per copper site. 



I. INTRODUCTION 

Our calculations were motivated in part by NMR ex- 
periments searching the spatial inhomogeneity of charge 
and spin density distributions in underdoped cuprates 
101 . These investigations are connected with the very 
important problem of high-temperature superconductiv- 
ity because the nature of a pseudogap is not completely 
understood. Therefore we continue our previous exam- 
ination of the singlet correlated band model for layered 
cuprates |^ with respect to its ability to describe the 
observed doping dependence of a pseudogap formation 
temperature [HJJ] and Fermi surface evolution 



II. BASIC RELATIONS FOR ONE BAND 
APPROACH 

We start from the Hamiltonian: 

+ X! dij^i^j + HcDW (1) 

where ^ ■q,'^,pd quasiparticle Hubbard- like oper- 
ators for the copper-oxygen singlet band, Jij is the su- 
perexchange constant of the copper spin coupling and 
Qij is a screened Coulomb repulsion of the doped holes, 
^ + ^ J2a '^T" + 2^-^ The quasiparticle inter- 
action HcDW mediated by the phonon field leads to a 
CDW transition ^. 

In addition to the usuall mean field approach we have 
taken into account that anticommutators of Hubbard-like 
operators can be affected by the doping index level per 
one unit cell, spin magnetization and charge modulation 
(non- Fermi statistics effect (NFS)): 

pr = + = ^ + (-i)^-'^^^ (2) 



The appearence of SDW we describe via the Fourier 
component: 

" ^ exp(ig^i?,), (3) 

where Qs is the instability wave vector with respect 
to SDW formation. CDW can be taken into account in 
the same wayr by introducing the Fourier component of 
doped holes: 

^1-^ exp(igi?,) (4) 

where qe is the instability wave vector with respect to 
CDW formation. It is easy to show that the homoge- 
neous part (So) of the doped hole function per one unit 
cell {6i) does not contribute to the expectation value Cq 
{i.e. Si — 5q + Ci). Generally the wave vectors of CDW 
(qe) and SDW (g^) can be different. Below we consider 
both as commensurate wave vectors qs = qe ^ (tt, tt) and 
incommensurate wave vectors qs = qe = (tt ± e^,, tt ± Ej,). 

Using Roth's decoupling procedure in the frame- 
work of a linear approximation we can write: 

< PfP- + vI>--^M.J'- > ^ 1+^ 4<g,., > _ 
< P/ > 2 ^ + (1 + 

-« I >+(-!)*-" <<»[l-i^l (5, 

where the angular brackets correspond to the thermo- 
dynamic expectation values. Then the equations of mo- 
tion are written as: 

^^^,Kf,-^^+,Kf,+vlK'' (6) 

where 
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<|-(-l)^-^5^>+G^+, + Gf (7) 
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and Efc is the energy dispersion in the normal state. 
The CDW order parameter is determined by the rela- 
tions il: 



k+q 



{l + 6o)N 



'9k' 



are dispersions of the lower Hubbard copper band and 
the singlet-correlated oxygen band, respectively and 



E 



dp 



E 



l + So 



'-k<nf+^,n^'>}^ (8) Ef = j2[ 



(1 - <5o) 



1-So 



(1 + So) 



< SiSj >]t}f e-Kp{ikRij) (16) 



< SiSj >]tlf exp{ikRi.j) (17) 



Gf =5^[A(a-,)-i3(c.,)x 



where 



{hujqf 9 {hup - |£fc|) 6 {hup - \ek - £fc+<?|) j 
{Sk - Sk+qf - {fl^qf 

Oq = AOmeV is the frequency of active phonon 



mode in the CDW transition. 

The thermodynamic values of the Fourier component 
< Cq > and < Sq > are calculated self-consistently: 

- ^ E[< > + < >] 

k 

< - 1 E[< > - < >] (10) 



describe their hybridization. If Sq goes to zero the spin 
dependent factors in E'^'^ and E^^ become zero due to 
(9) antifcrromagnetic correlations. At the same time, as one 
can see from the expressions for E^"^ and i?^^ the antifcr- 
romagnetic correlations enhance the interband coupling. 
So, at small doping level the energy dispersion will be 
determined by the expression 

ek{S^O) = '-^ + 

+ \^{ep- Sd)'^ + 16[ii2(cos fc, + cos ky)f (18) 



The correlation function is determined by 

"^k+q 



<*r4*r' >=p 



^Ik - E2k 



f{Elk)-f{E-2k)] (11) 



with 

Eik,2k 



'-^^^^±l[iSk-ek + qr+i\vUq 

(12) 



where is a transfer integral between nearest neigh- 
bor copper sites. From (18) one can obtain that the bot- 
tom of the band corresponds to the point (§, f ) of the 
Brillouin zone. This result coincides with experimental 
data We stress here that in contrast to Chubukov's 
model of the normal state of lightly doped cuprates ||] 
or slave boson approximations [Til ] long range order 
is unnecessary in our two band approximation. The ex- 
pression (18) is valid at < sf >= 0, i.e. for the usual 
paramagnetic phase. 



III. ENERGY DISPERSION AT LOW DOPING 



The one band approach is usually valid for a large 
enough doping level. Therefore if 6 goes to zero we ex- 
plored the two band model for layered cuprates |^,^,|lO|. 
In this case the energy dispersion is: 



£k 



+ EIP 1, 



-fc ^-fc , ^[^Ef - E^P)^ +AE^^PEf]^^^ (13) 



where 



1 (5 2 

Ek'^ £d + Et^r^ + 71 ~ < exp(iA:i?y, 



2 (1 - So) 



(14) 



^r=^p+E[ 



l + So 



(1 + So) 



< SiSj >]tij eicp{ikRij) 
(15) 



IV. NUMERICAL RESULTS AND DISCUSSION 



The spectral weight of the singlet band changes with 
doping level S (hole concentration per one unit cell of bi- 
layer) as 26/ (1 -I- 6). Thus, the condition of half-filling 
(which approximately corresponds to the optimal dop- 
ing level ) yields Sopt = 1/3. Because the bilayer unit cell 
contains two copper sites we conclude that Tc has a max- 
imum near hole concentrations a; « 1/6 per one copper 
site. 

The evolution of the Fermi surface for three doping lev- 
els are presented in Fig.l. With decreasing doping level 
the Fermi surface shrinks to the pockets near the point 
(■|, -1) of the Brillouin zone. 
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FIG. 1. Evolution of the Fermi surface of the singlet 
band for three doping levels; x — 0.025 (solid curve), 
X = 0.09 (dash-dotted) and x = 0.17 ( dashed ). 

This behavior qualitatively agrees with photocmission 
experimental data ||,^ and provides a good basis for 
the analysis of Peierls like instabilities, which as it was 
pointed out by many authors [|[|2|jl^, is very sensitive 
to the topological properties of the Fermi surface. 

Numerical solutions of the equations (8)- (10) are 
shown in Fig. 2. 
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FIG. 2. Critical temperatures vs. doping; (CDW )- 
long dashed curve, (SDW)- solid curve, T* (CDW)- 
dashed-dotted, (SDW)- short dashed. The parabola 
corresponds to the superconducting transition tempera- 
ture Tc (schematically). 

The values of hopping integrals in the one band ap- 
proximation were choosen between first, second and third 
neighbors as ti = 72.5, t2 = and — 12 (in meV) re- 
spectively. The set of other used parameters is presented 
in Table 1. 



Table 1. The values of spin-spin correlation function and 
superexchange plus Coulomb screening integral at different 
doping level 



X 


< S^Sj > 


Jo + Go (mcV) 


0.20 


-0.0381 


160 


0.15 


-0.0728 


225 


0.10 


-0.1287 


320 


0.05 


-0.2301 


320 



The short range interactions (superexchange and 
screened Coulomb repulsion), as a rule, yield the d- 
wave order parameters of CDW and SDW, whereas NFL- 
effects, together with phonon mediated interaction, lead 
to a anisotropic s-wave pseudogap component. In general 
the calculated CDW and SDW order parameters have 
s + id symmetry and the temperature dependence of s- 
and d-components are different. The critical temperature 
of the d-component (TJ) is always higher than T* at all 
considered doping levels. Both types of solutions (CDW 
or SDW) display the correct doping dependence. In ac- 
cording to photocmission data and NMR Q the 
critical temperature of the pseudogap goes down when 
the doping decreases. The phase diagram of supercon- 
ductivity (parabolic line) is given schematically as in . 
From Eqs. (8), (11) it is clear that the critical tempera- 
ture of the d- component CDW (TJ) is not sensitive to 
the external magnetic field. This result agrees well with 
recent experimental observation jist . Calculated mean 
field critical temperatures of CDW (TJ) are higher than 
for SDW in the complete doping range. This result is also 
consistent with the widely accepted opinion that a tran- 
sition towards the so called " stripe" phase in underdoped 
cuprates are charge rather than spin driven pq]. 



In Fig. 3 we show the magnitude of the instability wave 
vector Q = (tt ± e^;, tt ± Sy) vs. doping index. We have 
found the maximum of Sy around xcdw 0.09 for CDW 
and xsDW ~ 0.11 for SDW. It would be interesting to 
check our theoretical conclusion experimentally. 



It is easy to see from Eq. (10) that for a d-wave 
order parameter with pure commensurate wave vector 
q = Q — (tt, tt) both expectation values < Sq > and 
< Cq > vanish, but they are finite when the wave vector Q 
becomes incommensurate. Because in general the order 
parameters are complex, the expectation values < Sq > 
and < Cq > have real and imaginary components. Its 
values at the superconducting transition temperature are 
presented in Fig. 4. 
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FIG. 3. The variation of incommensurability compo- 
nent ey vs. doping index; for CDW- solid and for SDW - 
dashed curve, respectively. 

The real component < > is about 0.1 — 0.15 at 
temperatures T >Tc and disappears at T*. 
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FIG. 4. Calculated values < Cq > and < Sq > for tem- 
peratures around T^; Re < e, >(CDW)-dashcd curve, 
Re < eq > (SDW) - dotted-dashed, Im < sq > (SDW) 
- solid, Im < Cq >(CDW) -long dashed, Im < Cq >(SDW) 
-dotted. 

The imaginary parts of < > and < > in gen- 
eral are about 0.01 — 0.05 and stay up to TJ" and then 
drastically vanish. 

V. CONCLUSION 

To conclude we have calculated Fourier amplitudes 
of CDW and SDW order parameters in underdoped 
cuprates. Our calculation provides an explanation for the 



experimentally observed evolution of the Fermi surface 
and the doping dependence of the pseudogap formation 
temperature. The charge instability preforms the spin in- 
stability and the incommensurability has the maximum 
at a; ~ 0.1 holes per one copper site. We hope our present 
calculations will help in better understanding many of the 
strange features in the shape and linewidth of NMR in 
layered cuprates. 
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